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ABSTRACT 

We examine the Butcher-Oemler effect and its cluster richness dependence in the largest sample studied 
to date: 295 Abell clusters. We find a strong correlation between cluster richness and the fraction of 
blue galaxies, /s, at every redshift. The slope of the relation is similar for all richnesses, but at a 

given redshift, /s is systematically higher for poor clusters. This is the chief cause of scatter in the /b 
vs. z diagram: the spread caused by the richness dependence is comparable to the trend in Jb over a 
typical redshift baseline, so that conclusions drawn from smaller samples have varied widely. The two 
parameters, z, and a consistently defined projected galaxy number density. A, together account for all 
of the observed variation in /s within the measurement errors. The redshift evolution of Jb is real, and 
occurs at approximately the same rate for clusters of all richness classes. 

Subject headings: galaxies: clusters - evolution 



L INTRODUCTION 

The Butcher-Oemler (BO) effect provided some of the 
first evidence for the evolution of galaxies and clusters. 
Butcher & Oemler (1978, 1984) found an excess of blue 
galaxies in high redshift clusters in comparison with the 
typical early-type population observed in the central re- 
gion of local clusters (Dressier 1980). The BO effect has 
also been observed in more recent studies, which indicate 
an even stronger evolution of the fraction of blue galaxies 
in clusters (Rakos & Schombert 1995, Margoniner & de 
Carvalho 2000). 

The BO effect is an indicator of evolution in the cluster 
galaxy population, which may result from changes in the 
morphology and star-formation rates of member galaxies 
with redshift. The fact that blue galaxies are more com- 
monly found at higher redshifts, and the observation of 
an apparent excess of SO galaxies in low redshift clusters, 
lead to the suggestion by Larson et al. (1980) of an evo- 
lutionary connection between SO and spiral galaxies. This 
idea can explain the population of blue galaxies observed 
by Butcher & Oemler as spiral galaxies seen just before 
running out of gas, and the disappearance of this popula- 
tion in more evolved, low redshift rich clusters (Dressier 
et al. 1997, Couch et al. 1998). In the last decade, high- 
resolution Hubble Space Telescope images have allowed the 
determination of the morphology of these high redshift 
blue galaxies. Dressier et al. (1994), Couch et al. (1994, 
1998), and Oemler et al. (1997) found that the Butcher 
& Oemler galaxies are predominantly normal late-type 
spirals, and that dynamical interactions and mergers be- 
tween galaxies may be a secondary process responsible for 
the enhanced star formation. A recent study by Fasano 
et al. (2000) indicates that as the redshift decreases, the 
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SO cluster population increases while the number of spi- 
ral galaxies decreases, supporting Larson's original idea of 
spirals evolving into SOs. 

A second factor affecting the number of blue galaxies 
may be their environment. The fraction of blue galaxies 
seems to be correlated with local galaxy density, and with 
the degree of substructure in the cluster. Many studies 
have shown that the blue galaxies lie preferentially in the 
outer, lower density cluster regions (Butcher & Oemler 
1984, Rakos et al. 1997, Margoniner & de Carvalho 2000). 
Also, Small et al. (1998) studied 10 massive clusters and 
found that the fraction of blue galaxies in these clusters is 
smaller than observed in Butcher & Oemler (1984) clusters 
at the same redshift range. However, an opposite trend is 
found in non-relaxed clusters with a high degree of sub- 
structure, where the fraction of blue galaxies is higher than 
that observed in regular clusters (Caldwell & Rose 1997, 
Metevier et al. 2000). 

The evolution of these galaxies is probably correlated 
with its environment in the sense that spirals might con- 
sume and/or lose their gas, evolving to SOs, while falling 
into the higher density cluster regions. Although the exact 
mechanism responsible for the BO effect is not completely 
understood (Kodama & Bower 2000), most authors seem 
to agree that the effect is real. A different idea is pre- 
sented by Andreon & Ettori (1999) who show that the 
mean X-ray luminosity and richness of Butcher & Oem- 
ler (1984) clusters increases with redshift, and argue that 
the increase in the fraction of blue galaxies with redshift 
may not represent the evolution of a single class of objects. 

We present observations of the BO effect in a large sam- 
ple of 295 AbeU (1958, 1989) clusters of aU richnesses, 
with no further selection on the basis of morphology or 
X-ray luminosities. This is important because all previous 
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studies were based on samples biased toward richer clus- 
ters, with some samples being further selected according 
to morphology and/or X-ray luminosities, so that a com- 
bination of different selection effects could be mimicking 
the observed fsiz) relation. Although our sample inher- 
its biases existent in the Abell catalog, it should contain 
a more representative variety of clusters (in terms of de- 
gree of substructure, richness, and mass) at each redshift 
than any previous sample, and because of its large size 
should allow a better determination of the relation. Any 
results ^driven by selection effects present in previous sam- 
ples might become apparent when compared with this one. 
We describe the data in Section 2, the BO effect analysis 
in Section 3, and our conclusions in Section 4. 

2. SAMPLE SELECTION AND INPUT DATA 

The data presented in this paper were obtained to cali- 
brate the DPOSS-II (the Digitized Second Palomar Obser- 
vatory Sky Survey, Djorgovski et al. 1999). It comprises 44 
Abell clusters imaged at the 0.9m telescope at the Cerro 
Tololo Interamerican Observatory (CTIO) (Margonincr & 
de Carvalho 2000, hereafter MdCOO), and 431 clusters ob- 
served at the Palomar Observatory 1.5m telescope (Gal et 
al. 2000, hereafter GOO, and in preparation). The CCD 
images were taken in the g, r and i filters of the Thuan & 
Gunn (1976) photometric system, with typical l-o" mag- 
nitude errors of 0.12 in g, 0.10 in r, and 0.16 in i at 
r = 20.0™. We have also observed 22 control fields in or- 
der to assess the background contribution. More details 
concerning the data reduction, photometry, and catalog 
construction can be found in MdCOO and GOO. 

From this original sample we excluded 120 clusters ob- 
served with a small field of view CCD, and 31 very low 
redshift {z < 0.05) clusters, for which only the core region 
can be observed with our 13' x 13' images. Also, 26 clus- 
ter fields which exhibited galaxy counts comparable to the 
mean background, and one cluster with a bright star occu- 
pying ~ 25% of the CCD region, were excluded from the 
analysis. Since 2 Abell clusters had repeated observations 
from GOO and MdCOO, our final sample comprises 295 clus- 
ters. It is important to note that at variance with previous 
studies, the sample presented in this work is representative 
of all richness class clusters (21% i? = 0, 50% R=l, 21% 
R = 2, and 8% i? > 3). No furthca- sek;ction regarding rich- 
ness, morphology or degree of sub-clustering was applied 
when determining our sample. Spectroscopic redshifts for 
77 clusters were obtained from the literature, and for the 
remainder photometric redshifts were estimated with an 
rms accuracy of approximately 0.04 using the methodol- 
ogy described in MdCOO. 

3. ANALYSIS OF THE BUTCHER-OEMLER EFFECT 

The BO effect can be measured by comparing the frac- 
tion of blue galaxies in clusters at different redshifts. The 
(g—r) vs. r color-magnitude (CM) relation was determined 
by fitting a linear relation to the red galaxy envelope using 
the same prescription as in MdCOO. Red envelopes were 
subjectively classified as well-defined or not by visual in- 
spection. Galaxies were defined as blue if they had {g — r) 
colors 0.2"* below the linear locus in the CM relation. We 
used the spectral energy distribution of a typical elliptical 
galaxy (Coleman et al. 1980) to derive k{z) corrections for 



our sample, and also corrected the data for extinction us- 
ing the maps of Schlcgel et al. (1998). Because; the; main 
purpose of this work is to study the evolution of galaxy 
populations, fainter galaxies play a crucial role. Whereas 
Butcher & Oemler were limited to brighter galaxies by 
their photographic data, our CCD sample allows us to 
probe significantly deeper, selecting galaxies with r mag- 
nitude between M* — 1 and M* + 2, inside a region of 
radius 0.7 Mpc around the cluster. This fixed linear size 
was chosen because our CCD images cover a field of radius 
~ 0.5 — 4.0 Mpc at z = 0.03 — 0.38, and we want to study 
the same physical region and the same galaxy luminos- 
ity range for the entire sample. We assume a cosmology 
with Ho = 67 Km s~^ Mpc~^, and Qo = 0.1, in which 
M; = -20.16 (Lin et al. 1996). 

Most (61%) of the clusters in our sample are at z = 
0.1 — 0.2, so that the 0.7 Mpc central region, and the entire 
luminosity range between M* — 1 and M* + 2 can be ob- 
served. Clusters at higher redshifts are limited at brighter 
absolute magnitudes and a correction was applied in order 
to compare the fraction of blue galaxies in these clusters 
with the rest of our sample. The lower redshift clusters 
suffer from the opposite problem, since the brighter galax- 
ies will be saturated, and we are also limited by the field 
of view. Details about these corrections can be foimd in 
MdCOO, the only difference being that in the present work 
we adopt a more conservative brighter limiting magnitude 
of M* + 2 instead of M* + 3. 

The blue and total counts contain a mix of cluster 
members and background galaxies. Background variations 
place a fundamental limit on the accuracy of Jb measure- 
ments: no matter how well one determines the mean back- 
ground from a number of control fields, the background 
estimate for any particular cluster is no more accurate 
than the scatter in the control fields. For each cluster, we 
measured the effect of this scatter by computing /b us- 
ing background corrections from the 22 individual control 
fields scaled to the area of the CCD and the CM relation 
appropriate to that cluster. The final /b is given by the 
median and the rms is based on the central two quartiles of 
the distribution. Simple propagation of errors through the 
equation /s = I"'""'-'""^""' """" would give an 

^total , cluster 'f^ total , back ground 

overestimate because nuue, background and fltotal, background 

are correlated. The usual assumption that the error in 
Jb is mainly due to Poisson statistics in the background- 
corrected cluster counts does not properly take into ac- 
count the background variations and errors thus derived 
are about three times smaller than ours when applied to 
the same data. 

The final fractions of blue galaxies for the entire sam- 
ple are shown in the upper panel of Figure 1. For those 
clusters with multiple observations, we have used the ob- 
servation with smaller a/g. Clusters with spectroscopic 
and photometric redshift measurements are indicated by 
solid and open circles respectively. The individual error 
bars are not presented to avoid confusion in the plot, but 
the median af^ is 0.071. The lower panel presents only 
clusters that have: (1) cr/g < 0.071, (2) spectroscopic red- 
shift measurements, and (3) a well-defined CM relation. 
The 1-afg errors are indicated for these clusters. In both 
panels, the solid lines indicates a linear fit derived from 
the clusters with z < 0.25 (/s = (1.24 ± 0.07)^ - 0.01, 
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with an rms scatter of 0.1 for the entire sample shown in 
the upper panel, and fg = (1.34 ± 0.11)z — 0.03, with 
an rms of 0.07 for the sample shown in the lower panel). 
The fits were done with the GaussFit program (Jefferys 
et al. 1988) taking into account the measurement errors 
in /b- The dashed line indicates the rms scatter around 
the fit, and while the upper panel shows a larger scatter, 
the two derived relations are the same within the errors. 
A clear trend of strong evolution with redshift is seen. A 
Xv test applied to the fit presented in the upper panel re- 
sults in 1.36 if only /b errors are considered, and 1.13 
when the uncertainties in the photometric redshifts arc 
also taken into account. Such small Xv is due however to 
large measurement errors in this sample. A higher Xi^ of 
1.72 is obtained for the subsample of clusters with smaller 
afg and spectroscopic redshift (lower panel of Figure 1). 

Richness is a natural second parameter which might 
cause the range in fs at given redshift. The left upper 
panel of Figure 2 shows the fB{z) diagram with symbol 
sizes scaled by iV, the number of galaxies between M* — 1 
and M* + 2, inside 0.7 Mpc, after background correction. 
Only clusters from the lower panel in Figure 1 and at red- 
shifts between 0.1 and 0.2, where no corrections needed 
to be applied to compute /s, are presented in this figure. 
The solid line is the best- fit linear /b(z) relation for the 
26 clusters shown, and the dashed lines indicate the rela- 
tions obtained when the sample is subdivided in two 13 
cluster samples according to richness. The slopes of the 
three relations are the same within the errors (0.86 ± 0.23 
for the entire sample, 0.96 ± 0.26 for the subsample of 
richest clusters, and 0.86 ± 0.56 for the poorest ones). Al- 
though the uncertainties are large, the rate of evolution 
is approximately constant for all richnesses, and there is 
clearly an effect, with richer clusters tending to lower Jb- 
The suggestion by Andreon & Ettori (1999) that the ob- 
served increase of /b with redshift is caused by missing 
poor clusters at higher redshift is therefore no longer ten- 
able, because any such selection effects would serve to de- 
crease the redshift evolution of /b. The evolution in /b is 
real, and takes place in all richness classes. 

To ensure that richness and redshift are not correlated 
in this subsample, we plot N vs. z in the right upper panel 
of Figure 2. The best-linear fit (TV = -(38.2 ± 106.2)z -h 
(89.5 ± 14.5)) is indicated by the solid line, and the dashed 
lines represent adding/subtracting \-a uncertainties to its 
coefficients. Richness and redshift are uncorrelated in this 
sample. 

To gauge the richness importance, we investigate 
/B(z,iV) relations with various richness dependences 
(iV-i, 7V-''/2, 7V-^ and TV"'^), and found that cor- 
relates slightly better with /b. The final best-fit relation 
for the data, when both redshift and richness are consid- 
ered, and taking into account the errors in /b and A/', is 
/b = (1.03 ± 0.25)z (388.3 ± 111.4)iV-2 - 0.04, with 
Xv = 0.99. The F test indicates with > 99% confidence 
{F^^ = 23.8, for 25 to 24 degrees of freedom) that rich- 
ness is responsible for the improvement observed in Xv 
Richness is therefore extremely important in determining 
/b and it is in fact enough to account for the scatter ob- 
served in a simple /b(-z) relation. In the lower left panel 
of Figure 2 we plot /b vs. z + 378.6iV^2, which represent 
an edge-on view of the best-fit plane /b(-z, N), and in the 



lower right panel we present the redshift dependence of a 
richness-corrected-/B: the redshift cvohition is clear. 

Although all 26 clusters shown in Figure 2 were used to 
compute the relations indicated by solid lines in the figure, 
the two clusters represented by open circles (Abell 520, 
and Abell 1081) seem to deviate from the trend. These 
clusters also have the largest ajsi (derived from the scatter 
in the 22 control regions), and if excluded from the fitting, 
the relation changes to /b = I-IS^ + 480.56A/'-2 - 0.06 
{Xu = 0.95), indicating a slightly stronger redshift depen- 
dence. 

Finally, we used 11 clusters with ROSAT X-ray lumi- 
nosities to check for correlations with /b, but this small 
sample did not allow any conclusions. There is a possible 
trend of increasing Lx with redshift in the sense found by 
Andreon & Ettori (1999) in BO's original sample. How- 
ever, there is no clear trend of /b with Lx , arguing against 
any contamination of the /b(-2) relation by selection ef- 
fects. A comparison of the X-ray luminosities of these 
clusters with their /b could provide new insight to the 
BO effect, and we are in the process of obtaining X-ray 
fluxes and upper limits from the RASS (Rosat All Sky 
Survey) for all of these clusters, as well as a larger sample 
of new clusters generated from DPOSS-II. 

We stress that these results are not to be directly com- 
pared with fractions of blue galaxies as originally defined 
in B084. The most important reason for differences is the 
magnitude range used to calculate /b . The number of blue 
galaxies increases at fainter magnitudes and it is therefore 
natural that our /b measurements are in general higher. 
Also, the usage of a rigid physical scale for all clusters, 
instead of regions determined individually for each cluster 
according to its density profile, should yield different /b 
estimates. When /b is computed using the same absolute 
magnitude range used by B084, we find signs of evolution 
that are stronger than originally suggested in their work, 
and which are consistent with recent works by Rakos & 
Schombert (1995) and MdCOO. Limiting oiu' sample at 
brighter absolute magnitudes however increases the cr/^ 
(and (Tjv) because of the smaller galaxy count statistics, 
and probes a smaller fraction of the cluster population. 

4. SUMMARY AND CONCLUSIONS 

We compute /b for 295 randomly selected Abell (1989) 
clusters, including galaxies as faint as M* -|- 2 to sam- 
ple a larger range of the luminosity function and provide 
better statistics in each cluster. The resulting /b shows 
a stronger trend with redshift than did /b as originally 
defined by B084, consistent with the idea that the BO ef- 
fect is stronger among the late-type spirals and irregulars 
which donnnatc; the galaxy populations at intermediate 
and lower luminosities. 

This is the first sample large enough to allow the study of 
the /b variations at a given redshift. A Xu test applied to 
a simple linear fB{z) relation shows that the scatter in /b 
at a given redshift is consistently larger than the measure- 
ment error, indicating a real cluster-to-cluster variation. 
We investigate the richness dependence of the Butcher- 
Oemler effect, and find a strong correlation between /b 
and galaxy counts in the sense that poorer clusters tend 
to have larger /b than richer clusters at the same redshift. 
The inclusion of poor clusters tends therefore to increase 
the slope of the /b {z) relation, and is another factor (to- 
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gather with the inclusion of fainter galaxies) responsible 
for the stronger evolution observed in this sample when 
compared to previous works based mostly on rich clusters. 

We show that the fraction of blue galaxies in a cluster 
can be completely determined by its redshift and richness 
(within the measurement errors). The evolution in /b with 
redshift is real, and occurs at approximately the same rate 
for clusters of all richnesses. 
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Fig. 1. — Fraction of blue galaxies in the magnitude range between M* — 1 and M* + 2, and within 0.7 Mpc from the 
center of the cluster. The entire 295-cluster sample is shown in the upper panel, and only clusters with (1) aj^ < 0.071, 
(2) spectroscopic redshift measurements, and (3) well-defined CM relation are presented in the lower panel. In each panel 
the solid line indicates a linear /b(z) fit derived from the clusters with z < 0.25, and the dashed line represents the rms 
scatter of the clusters around the fit. 
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Fig. 2. Subsamplc of clusters with 0.1 < z < 0.2 from the lower panel in Figure 1. The solid line in each panel represents 
the best-fit to the entire sample. Left upper panel: ./b(z) diagram with marker sizes scaled by N (galaxy counts). The 
dashed lines indicate the fits computed separately for the richest (circles) and poorest (squares) clusters. Right upper 
panel: N vs. z. The linear N{z) fit is indicated by the sofid line, and the dashed lines represent adding/subtracting 
1-C7 uncertainties to its coefficients. Lower left panel: /b vs. z + 378.6N~'^ (edge-on view of the best-fit plane Jb = 
(1.03 ± 0.25)2;+ (388.3 ± 111.4)A''~^ — 0.04 to the data). Lower right panel: richness-corrected Jb as a function of redshift. 



